General partonic structure for hadronic spin asymmetries 

M. Anselmino, 1 M. Boglione, 1 U. D'Alesio, 2 E. Leader, 3 S. Melis, 2 and F. Murgia 2 

1 Dipartimento di Fisica Teorica, Universitd di Torino and 
INFN, Sezione di Torino, Via P. Giuria 1, 1-10125 Torino, Italy 
2 INFN, Sezione di Cagliari and Dipartimento di Fisica, Universitd di Cagliari, 
CP. 170, 1-09042 Monserrato (CA), Italy 
3 Imperial College London, Prince Consort Road, London SW7 2BW, U.K. 

O ' 

(SI ! (Dated: February 2, 2008) 

§ ! Abstract 
t-^ . 

\Q ' The high energy and large px inclusive polarized process, (A, Sa) + (B, Sb) — ► C + X, is con- 

sidered under the assumption of a generalized QCD factorization scheme. For the first time all 

<n ; 

transverse motions, of partons in hadrons and of hadrons in fragmenting partons, are explicitly 

m 
o 

On 
O 
IT) 

Several new spin and k± dependent soft functions appear and contribute to the cross sections and 



a: 



taken into account; the elementary interactions are computed at leading order with noncollinear 
exact kinematics, which introduces many phases in the expressions of their helicity amplitudes. 



to spin asymmetries; we put emphasis on their partonic interpretation, in terms of quark and gluon 



polarizations inside polarized hadrons. Connections with other notations and further information 
are given in some Appendices. The formal expressions for single and double spin asymmetries 
are derived. The transverse single spin asymmetry ^4^, for p* p — ► 7r X processes is considered 

X. 

in more detail, and all contributions are evaluated numerically by saturating unknown functions 
with their upper positivity bounds. It is shown that the integration of the phases arising from 
the noncollinear kinematics strongly suppresses most contributions to the single spin asymmetry, 
leaving at work predominantly the Sivers effect and, to a lesser extent, the Collins mechanism. 

PACS numbers: 13.88.+e, 13.60.-r, 13.15.+g, 13.85.Ni 
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I. INTRODUCTION AND FORMALISM 



There is, at present, no completely rigorous theory of single spin asymmetries in hadron- 
hadron collisions and inclusive particle production. Rigorous results about how different 
physical processes, including hadronic ones, are related to each other via factorization, only 
exist for the restricted case of collinear kinematics. But it is precisely in this kinematic 
situation that one cannot generate single spin asymmetries at leading twist. Thus, the 
introduction of intrinsic k± is crucial for a model of single spin asymmetries and we are 
therefore forced to rely on an intuitively reasonable calculational approach, within QCD, 
assuming a simple factorization scheme. This effectively neglects the role of the soft factors 
related to the Wilson lines which occur in the rigorous definition of k± dependent parton 
densities and fragmentation functions. 

In recent papers [1, 2] we have discussed such a formalism to compute cross sections for 
polarized and unpolarized inclusive processes, AB — > C X, fully taking into account parton 
intrinsic motion in distribution and fragmentation functions, as well as in the elementary 
dynamics. In particular, in Ref. [2] the emphasis was on the importance of the many phases 
appearing in the computation of helicity amplitudes in noncollinear configurations, and their 
role in suppressing the contribution of the Collins mechanism [3] to transverse single spin 
asymmetries. Many other contributions to polarized and unpolarized cross sections, and to 
single and double spin asymmetries, were not discussed, referring to a later paper for the 
full treatment of the most complete case. 

We consider here such a general case. Let us start from Eq. (8) of Ref. [2]: 



E c d g( A ' s A)+(B,s B )^c + x ^ dx a dx b dz 
d 3 Pc abcd{x} J WXaX^S 



E / dx ; dXbdz 2 d 2 k ±a d"k ±b d*k ±c 5(k ±c ■ p c ) j(k ±c ) 



X Px^Jt fa/A,S A (X a , k ±a ) p^f* f b/B ,S B (x b , k ±b ) (1) 

111 u, o' o 

which gives the cross section for the polarized hadronic process (A, Sa) + (B, Sb) — > C + X 
as a (factorized) convolution of all possible hard elementary QCD processes, ab — > cd, with 
soft partonic polarized distribution and fragmentation functions. In Eq. (1) s, t and u are 
the Mandelstam variables for the partonic reactions and the detailed connection between 
the hadronic and the partonic kinematical variables is given in full in Appendix A. 

Let us clarify the physical meaning of Eq. (1) - our starting point - by making detailed 



comments on its notation and contents. 



• A and B are initial spin 1/2 hadrons (typically, two protons), in pure spin states 
denoted by Sa and Sb respectively, with corresponding polarization vectors P A and 
P B (notice that P A ' B are actually pseudovectors) . We set Sa,b = for unpolarized 
hadrons (P A ' B = 0). Ec and p c are the energy and three-momentum of the final 
detected particle (typically, a pion). Throughout the paper, we work in the AB cm. 
frame, assuming that hadron A moves along the positive Z cm -axis and hadron C is 
produced in the (XZ) cm plane, with (pc)x cm > 0. We define as transverse polarization 
for hadrons A and B the F cm -direction, often using the notation 

t for Pd = 1 and PS = -\ [ for P$ = -1 and PB = 1 . (2) 

1 * cm -* cm ^ * cm -* cm V / 

The longitudinal spin states are labelled by their helicities: \a,b = ±1/2 (sometimes 
just written as ±) corresponding to Pz cm = ±1 and to P£ m = =Fl respectively. The 
opposite signs for hadrons A and B originate from the fact that their helicity frames, 
as reached from the overall cm. frame, have opposite Y and Z axes [4], see Eq. (D3). 
The general case of hadrons transversely polarized along a generic direction <ps A in the 
(XY) cm plane is treated in Appendix B. 

• The notation {A} implies a sum over all helicity indices. x a , Xb and z are the usual 
light-cone momentum fractions, of partons in hadrons (x a ^) and hadrons in partons 
(z). k± a (k±b) and k±c are respectively the transverse momenta of parton a (b) with 
respect to hadron A (B), and of hadron C with respect to parton c. We consider all 
partons as massless, neglecting heavy quark contributions. 

• With massless partons, the function J is given by [1] 

2 



Ec + VP& " k ±c 

• P\ A w A 1S the helicity density matrix of parton a inside the polarized hadron A, with 
spin state Sa- Similarly for parton b inside hadron B with spin Sb- Notice that 
the helicity density matrix describes the spin orientation of a particle in its helicity 
frame [4]; for a spin 1/2 particle, Tr (c^p) = Pj is the i-component of the polarization 



vector P in the helicity rest frame of the particle. Obviously, for a massless parton 
there is no rest frame and the helicity frame is defined as the standard frame [4] in 
which its four-momentum is = (p, 0, 0,p) (see also Appendix D). f a /A,s A ( x a, k± a ) is 
the distribution function of the unpolarized parton a inside the polarized hadron A. 
We shall also denote by fgjSj ^ ne number densities of partons a, with spin along the 
i-axis, inside a hadron A with spin along the J-axis: i — x,y,z stand for directions in 
the parton helicity frame, whereas J = X, Y, Z refer to the hadron helicity rest frame. 

• The MA c ,A d ;A a ,A 6 's are the helicity amplitudes for the elementary process ah — > cd, 
normalized so that the unpolarized cross section, for a collinear collision, is given by 

da ab -> cd 1 1 v ~ 2 

dt ~ 16^4 . f-f . 1 A o,A d ;A a ,A 6 | • ( ) 

"A A' 

• D x c ' x ,°(z, k±c) is the product of fragmentation amplitudes for the c ^ C + X process 

Dx^if = f^ XXx ^a c ,a x ;A c V* x , c:Xx . x , c , (5) 

where the Y] stands for a spin sum and phase space integration over all undetected 
Tx,x x 

particles, considered as a system X. The usual unpolarized fragmentation function 
Dc/ C {z), i-e. the number density of hadrons C resulting from the fragmentation of an 
unpolarized parton c and carrying a light-cone momentum fraction z, is given by 




Eq. (1) is written in a factorized form, separating the soft, long distance from the hard, 
short distance contributions. The hard part is computable in perturbative QCD (pQCD), 
while information on the soft one has to be extracted from other experiments or modeled. As 
already mentioned and discussed in Ref . [2] , such a factorization with noncollinear kinematics 
has never been formally proven. Indeed, studies of factorization [5-7], comparing semi- 
inclusive deep inelastic scattering (SIDIS) and Drell-Yan reactions have indicated unexpected 
modifications of simple factorization, and the situation for inclusive particle production in 
hadron-hadron collisions is not yet resolved. Thus, our approach can only be considered as 
a reasonable phenomenological model. Of course, the perturbative calculation of the hard 
part is only reliable if the hard scale - in this case the square of the transverse momentum 
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of the final hadron, p\ - is large enough; in our case -p\ > 2.25 (GeV/c) 2 . It turns out that 
the data demand [1] an average value of k\ ~ 0.64 (GeV/c) 2 for the intrinsic transverse 
momentum of the parton distributions. This is relatively small compared to 2.25 (GeV/c) 2 , 
but complications can arise from the tail of the Gaussian distribution, as was discussed in 
Ref. [1] and will be commented on in Section V. 

The intrinsic motion arises both from parton confinement and from QCD gluon emission. 
In that, our approach, based on perturbative computations performed at LO in the strong 
coupling constant, with noncollinear kinematics, could partially and effectively contain some 
of the effects related to soft gluon emissions and the threshold resummation of large log- 
arithmic perturbative corrections, recently performed within proper collinear factorization 
[8]. A study of weighted single spin asymmetries for double-inclusive production in hadron- 
hadron collisions, based on k± factorization using a diagrammatic approach, has appeared 
very recently [5]. 

In the next Section we discuss in detail the soft contributions to Eq. (1), related to 
parton distribution and fragmentation functions, while in Section III we give the explicit 
analytical expressions of all elementary amplitudes, convoluted with the corresponding soft 
functions. Some contributions to the unpolarized cross section and the transverse Single 
Spin Asymmetry (SSA) are analytically discussed in Section IV. Numerical estimates of 
the maximal contributions of the different spin mechanisms, both to the cross section and 
the transverse SSA, are presented and discussed in Section V. General conclusions and 
comments are given in Section VI. Finally, the full noncollinear partonic kinematics and its 
relation with the overall hadronic variables is discussed, for convenience and completeness, 
in Appendix A; the formal relationships between the hadron and the parton polarization 
are widely studied in Appendix B, and the connection with other formalims is explicitely 
worked out in Appendix C. Useful definitions of helicity frames are given in Appendix D. 

II. SOFT PHYSICS 

Although Eq. (1) has already a clear physical interpretation, we would like to express the 
parton density matrix elements in terms of parton polarizations, so that, when performing 
the helicity sums, each term has a direct partonic meaning. 

Notice that the parton polarizations are, of course, related to their parent hadron polar- 
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izations. The way the hadron spin is transferred to the partons can be formally described, 
in general, by bilinear combinations of the helicity amplitudes for the process A — > a + X 
(distribution amplitudes) [2, 9]. Therefore, one could equally well interpret Eq. (1) either 
in terms of parton polarizations or in terms of the distribution amplitudes. We follow here 
the former approach, which is somewhat more direct. However, the latter approach offers 
a deeper understanding of some of the basic properties of our factorized scheme (e.g. the 
parity properties) and allows a direct comparison with other formalims used to describe the 
same spin effects. In Appendix B we give the full correspondence between parton polariza- 
tions and the distribution amplitudes, and in Appendix C we derive the explicit relations 
between our formalism and that of the Amsterdam group [10]. 



A. Quark polarizations 



The helicity density matrix of quark a can be written in terms of the quark polarization 
vector components, P a = (P", P y a , P") = (P$ cos Sa , P$ sin (p Sa , P£) , as 



a/A,S A _ 



( 



>A,S A 



1 + P a P «- iP « 

P? + iP* 1 - P? 



( 



1 + Pf PZ e -< * 4 « 



'A,S A 



pa c i<frs 



1 - Pf 



(7) 



'A,S A 



where, as explained above, the x, y and ^-directions are those of the helicity frame of parton 
a. Eq. (7) satisfies the well known general properties: 

= 1 (8) 



p;+ + p-- 
p++ - p a - 



= p° = pf 



2 Rep a _ + = 2 Rep a + _ = P£ = P? cos <f) Sa 
2 Imp" = -2 Imp^_ = P„ a = P£ sin (f) Sa . 



(9) 
(10) 

(11) 



When performing the sum over the helicity indices A a , A^ and A 6 , X' b in Eq. (1), one obtains 
products of terms of the form 

(Pj fa/A,S A ) = Jsj/Sa ~ f-Sj/S A = ^fsj/SA ' (I 2 ) 

where j = x, y, z. Similarly for parton b inside hadron B. We use the notations: 

(pf f a/ A,s Y ) = m: ]/Sy = /;,/t - /-,,/t = A^ J/t (x a , k ±a ) (is) 



(Pi fa/A, 



/A,S Z , 



Ml 



Sj/Sz 



fsj/+ ~ f-sj+ = Afsi/+( x a, k 



Sj/+ 



Sj /+\ 



(fa/A,S Y ) = fa/A(x a , k± a ) + ^ ^fa/S Y ( x a, k±a) ■ 



(14) 
(15) 



These amount to eight independent quantities, which represent the (k± unintegrated) distri- 
bution functions of partons a (= q, q) with polarization P a (defined in the partonic helicity 
frame) inside hadron A with spin Sa (specified in the hadronic helicity frame). All of these 
functions have a simple direct physical meaning: for instance, the x-component of Eq. (13) - 
(Px fa/A,s Y ) ~ represents the amount of polarization along the x-axis (in the partonic helicity 
frame) carried by partons a inside a transversely polarized hadron (A,Sy)] (Py f a /A,s Y ) * s 
related to the k± dependent transversity distribution, which, upon integration over d 2 kj_, 
gives the familiar transversity function h\(x) or A T q(x) (see also Appendix B). Similarly, 
the z-component of Eq. (14) - (P" f a /A,s z ) ~ is the unintegrated helicity distribution, which, 
once integrated over the transverse momentum, gives the usual helicity distribution Aq(x) 
or g\{x). 

Notice that two independent distribution functions appear in the definition of f a /A,s Y i 
which is the only term in the sum over A a , which corresponds to parton a being un- 
polarized: f a /A(x a ,k± a ), the unintegrated number density of unpolarized partons a inside 
the unpolarized proton A, and Af a /s Y , the Sivers function [11]. The latter permits the 
number density of unpolarized partons a to depend upon the transverse polarization of the 
parent hadron A. In general, for a hadron A in a pure spin state Sa and corresponding unit 

* A 

polarization vector P , one has: 

Af a /s A (X a , k± a ) = fa/S A (x a , k ±a ) - f a /-S A (x a , k ±a ) 

= A N f a/A r (x a , k ±a ) (pa x k ±a ) ■ P A . (16) 

In the last term of the above expression we have explicitly extracted the angular dependences, 
according to the so-called "Trento conventions" [12]: p A is the unit vector along the hadron 
A three-momentum, k± a = \k± a \ and k± a = k± a /k± a . Parity invariance allows to have a 

" A 

non-zero Sivers function only for transverse spin, p A P =0. Often A N f a / A t {x a , k± a ) alone 
is referred to as the Sivers function (see Appendix B for related expressions) . For a generic 

"A 

transverse polarization direction P = (cos0s A , sin05 A , 0), one has Af a /s A {x a) k^ a ) = 
A N fa/At ( x a, k± a ) sin((j) s A — 4> a) , where (p a is the azimuthal angle (in the hadronic cm. frame) 
of k ±a . 

According to our configuration the hadron transverse polarization is chosen along the 
+F-direction (f); notice that Y = Y cm for the hadron moving along the +Z cm -direction, 
while Y = —Y cm for the hadron moving along — Z cm) as already noticed after Eq. (2). Then, 
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Eq. (16) reads: 

A/a/Sy (X a , k± a ) = fa/S Y (x a , k± a ) - f a /-S Y (x a , k± a ) 

= A N f a/A i(x a ,k ±a ) COS0 a . (17) 

Similarly, the Boer-Mulders mechanism [10, 13] (see Appendix B) allows partons to be 
transversely polarized inside an unpolarized parent hadron. In general, this can be expressed 
by: 

Pj fa/A = fsj/A( X a, k± a ) ~ f- Sj /A( X a, k± a ) = Af^/ A (x a , kj_ a ) 

= & N f a y A (x a) k ±a ) (p A x k ±a )j , (18) 

where P? is the j-component of the parton polarization in the parton helicity frame (j = 
x,y,z). The above equation can also be written as [12] 

Afs/A( X a, k± a ) = fsl A {Xai k± a ) — fl s / A (x a , k± a ) 

= A N f a y A (x a , k ±a ) (p A x k ±a ) ■ P a , (19) 

* a 

where s and P denote respectively a generic parton spin state and the corresponding unit 
polarization vector, in the parton helicity frame (as reached from the parent hadron helicity 
frame). Notice that, according to our configuration, in the hadronic cm. frame y points 
along the Z cm x k± a direction, Eq. (D4). It follows that for nucleons moving respectively 
along the ±Z cm -direction one has 

&f: y/A (Xa, fc±o) = ±* N fat/A ( x a, fc±a) ■ (20) 

It also follows that the analogous function for the x-direction is zero, Af^^ A = 0. The 
function A N f a y A (x a , k^_ a ) alone is often referred to as the Boer-Mulders function. 

Moreover, one can show that the Boer-Mulders function is the same which appears in the 
y-component of Eq. (14), 

Nl/A = (Py a fa/A) = fa/A,S z ) = A£ /+ , (21) 

due to parity invariance. 

It is worth mentioning that the function {Py f a /s Y ) = ^fs y /s Y = ^fsy/1 can ^ e decom- 
posed into two terms, the Boer-Mulders term which is independent of the hadron transverse 
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polarization, and a term which changes sign when the hadron polarization direction is re- 
versed: 

Afsy/Sy = ^fsy/A + A fs y /S Y > ( 22 ) 



with 



Notice that 



A fs v /S Y = o [A/s„/T ~ A/",,] - _A fsyl-Sy ■ ( 23 ) 



^ft/S Y - A /s"/Sy - ~Afs x /-S Y ■ ( 24 ) 

B. Gluon polarizations 

Let us now consider the gluon sector (a first study of the unintegrated gluon distribution 
functions can be found in Ref. [14]). The helicity density matrix for a massless particle with 
spin 1 can be written as 



g/A,S A = 



1 1 + P z g V - \ 1 



\ T ° + iT 2 1 P * ) A ,S A \ llinV 1 larc JA,S A 



( 1 _|_ p9 _p9 p -2i<t> 
~ circ 1 lin 

_p9 2i</> i _ pS 



(25) 



and we consider it for a gluon g inside the hadron A, in a spin state Sa- Eq. (25) refers, 
in general, to a mixture of circularly and linearly polarized states. P^ rc corresponds to Pf, 
the gluon longitudinal polarization. The off-diagonal elements of Eq. (25) are related to the 
linear polarization of the gluons in the (xy) plane at an angle <fi to the x-axis. The x, y and 
z axes refer to the standard gluon helicity frame, in which its momentum is p^ = (p, 0, 0,p). 
Pf in is expressed in terms of the parameters Tf and ■, which are closely related to the 
Stokes parameters used in classical optics; they play a role formally analogous to that of 
the x and y-components of the quark polarization vector in the quark sector. The use of 
the parameters Tf and makes the gluon distribution functions formally similar to those 
for the quarks and considerably simplifies all the formulae for the spin asymmetries given in 
Sections III and IV. 

In analogy to the quark helicity density matrix, Eq. (25) shows that: 

p 9 ++ + p s _ = 1 (26) 

p 9 ++ - f!_ = PS = P? irc (27) 

2 Rep 9 _ + = 2 Rep%_ = T° = -Pf m cos(20) (28) 

2 Imp 9 + = -2 \mp 9 + _ = 7? = -P? m sin(20) . (29) 



As for the quark sector, there are eight independent gluon distribution functions, which, 
following Eqs. (13-15), we label as 

(Vf g/ASY ) = Aft in (x g ,k ±g ) (30) 

(T 2 9 f g/ASY ) = Af^ n (x g ,k ±g ) (31) 

(PS W) = A/f z/Sy = fl n - f_ Szn ee Af szn (x g , k ±g ) (32) 

(Vf g/ASz ) = Af^ /+ (x g ,k ±g ) (33) 

(T 2 9 f g/ASz ) = Af^ /+ (x g ,k ±g ) (34) 

(P? f 9 /A,s z ) = Mlis z = fl/ + ~ ts z/+ = A/i /+ (x p , k ±g ) (35) 

(fg/A,S Y ) = fg/A(x g , k ±g ) + - Af g/ s Y (x g , k ±g ) . (36) 

Notice that Af^ + (x g ,k± g ) is the usual k± g dependent gluon helicity distribution func- 
tion Ag(x g , k± g ). The interpretation of Afc^/SA as difference of linearly polarized gluon 
distributions is discussed in the sequel and in Appendix B. 

In analogy to Eqs. (22) and (23) we also define a new quantity which changes sign when 
the hadron polarization direction is reversed [see Eq. (B39)]: 

Af% n =Af\ /A + A-f% /v (37) 

with 

A ^/T = \ [ A /r l/T " A/J/J = -A-/f i/A . (38) 
Although gluons cannot carry any transverse spin, there is a strong analogy between 
transversely polarized quarks and linearly polarized gluons; for example, analogous to the 
Boer-Mulders case for quarks, it is possible to have a linearly polarized gluon inside an 
unpolarized nucleon, corresponding to a non-vanishing Tf f g /A = Af^^ A . This mechanism 
has never been explored before. Its structure is linked to the spin 1 Cartesian tensor Ty (see, 
e.g., Section 3.1.12 of Ref. [4]), which is symmetric and traceless. For a massless particle 
one has: 

T. = Jg (39, 

%. = M{T* X -T m ) % = 2^T„ V . (40) 

Because of (39), the traceless condition and parity invariance, it is only possible to construct 
one scalar structure that depends non-trivially on the T^. 
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Using the three- vectors at our disposal - the gluon momentum p, its transverse momen- 
tum k± and the parent hadron momentum p A - we define: 

k±- (k ± -p)p „ 
w = - v ^ f^i: v = p x fc ± (41) 

PPa 

and introduce a tensor T whose components are T^. The only possible structure is then: 

~ "I 



Tij f g/A (x,k ± ) = ^| ^ A N f°- i/A (x, k ± ) (uiUj - ViVj) 

(42) 



^ /s/a^, fc±) + - ZPiPj) 



which is the gluon tensorial analogue of Eq. (18). When nucleon A moves along or opposite 
the Z cm -axis this reduces to: 



Vf 9 /A{x, k ± ) = Afl /A (x, k ± ) = A N fl /A (x, k ± ) (43) 

in analogy to Eq. (20). Notice that, in this case, there is no ± sign on the r.h.s. of Eq. (43). 

One can also show that the linear polarization Tf is independent of any longitudinal 
polarization of the nucleon, i.e.: 

as in Eq. (21). 

C. Quark and gluon fragmentation functions into unpolarized hadrons 

As already mentioned in Section I, for the fragmentation process in general we define 

D%$>(z, k ±c ) =4 A c ,x x ;xM VI, x x . tX ,{z, k ±c ) . (45) 

The analogous quantity for parton distributions can be found in Eq. (B2). T>\ ,a x ;A c is the 
fragmentation amplitude describing the process c — > C + X in which the parton c from 
the elementary scattering ah — > col generates the detected final hadron C, with light-cone 
momentum fraction z and transverse momentum k± C - If we denote by 0^ the azimuthal 
angle of the hadron C in the parton c helicity frame, we have 

£a c ,a x; a>, k ±c ) = V Xc , Xx , Xc (z, k ±c ) e^o , (46) 
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similarly to Eq. (B4) for parton distribution amplitudes. Eqs. (45) and (46) then give the 
generalized fragmentation function 

£%$>(z, k ±c ) = D X £ x f(z, k ±c ) e^-^c , (47) 

while the corresponding generalized distribution function is given in Eq. (B5). 

If hadron C is unpolarized, the generalized fragmentation function D simply becomes 

D C ^ K (z, k ±c ) = £ D X X ± (z, k ±c ) = D C X ^ (z, k ±c ) e^-^c , (48) 

and fulfills the parity properties given by 

D%_ K (z, k ±c ) = (-1) 2 - (-l)V^ D C x ^,(z, k ±c ) ■ (49) 

If parton c is a quark, s c — 1/2 and the helicities A c and A^, will be either = +1/2 or —1/2, 
whereas if parton c is a gluon, s c = 1 and A c and A^, will be either = +1 or —1. 

Quark fragmentation functions 

For quarks, from Eqs. (45), (48) and (49) we obtain the following relations 

D+{ q (z, k LC ) = D% q (z, k ±c ) = D c/q {z, k ±c ) 

D C l q (z,k LC ) = D^(z,k ±c ) (50) 
for equal helicity indices, and 

D%\z, k ±c ) = D% q (z, k ±c ) e^c = -D% q (z, k ±c ) e«% 

D% q (z, k ±c ) = D% q (z, k ±c ) e-«% = -D c ; /q {z, k ±c ) (51) 

[D°L q (z,k ±c )}* = -D° / _ q (z,k ±c ) 

for unequal helicity indices. 

Dc/ q (z, k±c) is the k±c dependent fragmentation function describing the hadronization 
of an unpolarized quark q into an unpolarized hadron C . Notice that it does not actually 
depend on the direction of k±c, but only on its modulus. When integrated over the intrinsic 
transverse momentum, this function gives us the usual unpolarized fragmentation function 
D C /c(z), see Eq. (6), 

Dc /q {z) = \Y.j d 2 k ±c D c x ^(z,k ±c ). (52) 
12 



Eqs. (51) tell us that the fragmentation function D + _(z,k±c) is an independent purely 
imaginary quantity It is related to the Collins quark fragmentation function by the following 
expression: 

-2iD°L q (z,k ±c ) = 21m D°L q (z,k ±c ) = A N D c/q ,(z, k ±c ) , (53) 

and gives the difference between the number densities of unpolarized hadrons C resulting 
from the fragmentation of a quark q polarized along the +y-direction and a quark polarized 
along the — y-direction, in the quark helicity frame in which the fragmentation process occurs 
in the (xz) plane. In general one has, analogously to Eq. (16) for the Sivers function, 

AD C / q , s (z, k± c ) = D c / q , s (z, k± c ) - D c/q _ s (z, k± c ) 

= A N D c/q1 (z, k ±c ) (P q x k ±c ) ■ P q ■ (54) 

If P points along the ^-direction Eq. (54) reads, in analogy to Eq. (17): 

AD c / q ,s y (z, k ±c ) = D c / q , Sy (z, k ±c ) - Dc/q,-s v (z, k ±c ) 

= A N D c/q t(z,k ±c ) cos0g, (55) 

consistently with Eq. (48). The explicit expression of 0^ in terms of the overall hadronic 
variables, in the AB cm. frame, can be found in Eq. (45) of Ref. [2] and in Appendix A, 
Eq. (A28). 

Gluon fragmentation functions 

The gluon fragmentation functions with equal helicity indices obey the same parity rules 
(50) as the quark ones: 

D% 9 (z, k ±c ) = D% 9 (z, k ±c ) = D c/9 (z, k LC ) 

D C/ - 9 (z,k ±c )=DZ 9 (z,k ±c ); (56) 

however, as implied by Eq. (49), a different sign, with respect to the quark case (51), appears 
in the parity relations for the generalized gluon fragmentation functions with unequal helicity 
indices: 

D°L'(z, k ±c ) = D c + L 9 (z, k ±c ) e 2 ^" = D% 9 (z, k ±c ) e 2 ^ 

D% 9 (z, k ±c ) = D% 9 (z, k ±c ) e- 2 ^c = D c/9 {Zj k±c) e -«#g (5?) 

[D% 9 (z,k ±c )}* = D% 9 (z,k ±c ). 
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The above equations show that D + _(z,k±c) is an independent, real quantity. Notice that 
the gluon Collins fragmentation function cannot exist, since there is not such an object as 
a transversely spin polarized real gluon. However, similarly to what happens for the gluon 
parton distributions, the fragmentation function D+^(z, k±c) is related to the fragmentation 
process into a spinless hadron C of a linearly polarized gluon. In analogy to Eq. (53) we 
have 

2D% g (z,k ±c ) = 2ReD c l/_ 9 (z,k L c) = A N D c/Tf (z, k LC ) , (58) 

which gives the difference between the number densities of unpolarized hadrons C resulting 
from the fragmentation of a gluon linearly polarized along the ^-direction and a gluon lin- 
early polarized along the y-direction, in the gluon helicity frame in which the fragmentation 
process occurs in the (xz) plane. 

III. KERNELS 

As we can see from Eq. (1), the computation of the cross section corresponding to any po- 
larized hadronic process (A, Sa) + (B, Sb) — > C + X requires the evaluation and integration, 
for each elementary process a + b — > c + d, of the general kernel 

E(S A , S B r b - cd = Y,ptt f*/A,s A (x a , k ±a ) p^ A f f b/B ,s B (x b , k ±b ) (59) 

{A} 

X M\ AAA My \ .y y D X x C 'y C (z,k ±C ) ■ 

Whereas the hadronic process (A, Sa) + (B, Sb) — > C + X takes place, according to 
our choice, in the (XZ) cm plane, all the elementary processes involved, A(B) — > a(b) + X, 
ab — > cd and c — > C + X do not, since all parton and hadron momenta, p a , p b , p c have 
transverse components k± a , k± C - This "out of (XZ) cm plane" geometry induces in the 
fragmentation process the phase given in Eq. (48) and, in the distribution functions, the 
phase appearing in Eq. (B5). 

Analogously, the elementary QCD process ab — * cd, whose helicity amplitudes are well 
known in the ab center of mass frame, is not, in general, a planar process anymore when 
observed from the AB center of mass frame, the laboratory frame, where we are performing 
our computations. However, we can go from the actual p a p b — > p c Pd configuration, as seen 
in the laboratory frame, to the canonical one in which the ab — > cd process takes place 
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in the ab cm. frame and in the (XZ) cm plane, by performing one boost and appropriate 
rotations, as described in full detail in Ref. [2]. These transformations introduce some highly 
non trivial phases in the helicity amplitudes Ma c ,a ;A q ,a , which are the direct consequence of 
the complicated non planar kinematics. The relation between these amplitudes (which we 
need in our computations) and the usual, canonical amplitudes M°, defined in the partonic 
ab — > cd cm. frame, is the following [2] 

M\ x \ \ = M? x \ , e -H.KZa+hZb-Kic-* d U) e -i[{K-\)L-(K-* d %\ e KK-\Wc (qq) 

with £j, £j (j = a, b, c, d) and cj)" defined in Eqs. (35-42) of Ref. [2] and in Appendix A. The 
parity properties of the canonical cm. amplitudes M° are the usual ones: 

^-A c ,-A d ;-A a ,-A 6 = VaVbVcVd(—]-) Sa+Sb ~ Sc ~ Sd ( — l) < - Xa ^ X b^ ( ^ X ^ X d) M^^.^^ , (61) 

where rji is the intrinsic parity factor for particle i. For massless partons there are only 
three independent elementary amplitudes M° corresponding to the ab — > cd processes we 
are interested in; this allows us to adopt the following notation 





= M°e iipi 








M_ +; _+ 


= M 2 °e^ 2 






(62) 




= M 3 °e^ 3 , 








as 










M°, 


,+ = M° 2 




-,- - , 


(63) 



and the phases ip±, ip 2 and ip s are given by replacing in Eq. (60) the appropriate value for 
the helicities X i: i = a,b,c,d. Indeed, the + and — subscripts refer to (+1/2) and (—1/2) 
helicities for quarks, and to (+1) and (—1) helicities for gluons. 

All other amplitudes are obtained from Eqs. (60), (62) and (63), exploiting the parity 
properties (61); notice that the presence of the phases tpj implies that the parity relations 
for the amplitudes M are not as simple as those for the M°. From Lorentz and rotational 
invariance properties [4] one can obtain the following useful expressions relating the canonical 
amplitudes for processes which only differ by the exchange of the two initial partons, a <-> b, 
or of the two final partons, c^> d: 

M° x ' b Vx\ {0) = M° x ' fTA (f - 9) e -^- A d ) (64) 
Ki7£^ = ^a°:"a 6 (- - 0) e-^~ x , (65) 
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where the scattering angle 9 is defined in the canonical partonic cm. frame. To be precise, 
the above relationships hold up to an overall, helicity independent, phase; since only bilinear 
combinations of the amplitudes occurr in the expressions for physical observables, we fix such 
a phase to be +1. 

There are eight elementary contributions ab — > cd which we have to consider separately 



q a qb -> q c qd , 9a9b -> 9c9d , 

qg ^qg, gq-+ gq, 
qg -»■ gq, gq^qg, 
g a gb^qq, qq^g c gd, 



(66) 



where q can in general be either a quark or an antiquark. The subscripts a, b, c, d for 
quarks, when necessary, identify the flavour (only in processes where different flavours can be 
present); for gluons, these labels identify the corresponding hadron (a — > A, b — > B, c — > C). 
By performing the explicit sums in Eq. (59), we obtain the kernels for each of the elementary 
processes. Note that the new aspect of our calculation is the oppearance of the phases 
which is a reflection of the noncollinear kinematics. For convenience we also give explicit 
expressions for the combination of the partonic cm. amplitudes M? which are needed. 

1) q a qb -> q c qd processes 



£(^,S B )*"™« = -D c/c (z,k ±c ) f a/SA (xa,k ±a ) fb/s B (x b ,k ±b ) x 

M«|2 + |M0|2 + |M0|2) + pa pb (j^0|2 _ |^0|2 _ |^0|2^j 



(p; p, b + p; j*) cos(^ 3 - y? 2 ) - (p; p* - p; p*) smfo - y? 2 ) 



+ 2M 2 ° M 3 ° 



- 2 A ^C/cT (Z, k ±C ) fa/S A (Xa, fe±o) /&/S s fab, X 

{m° M 2 ° [P x a sin(^ - y? 2 + 0c) " Py cos(v^i " V?2 + 0c) 
+ Ml Ml \P b x sin( V 9 1 - y? 3 + 0c) - cos (^i -V3 + 0c) 



(67) 



where (including color factors) 



l M i°| 2 = 9^ 



+ 5, 



ab 



3 t-u , 



|M 2 °| 2 



9 9s i 2 
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1 2 



A 8 ^ 



sit 1 s 



M°M$ = 5 ab -gU^-~ 



st Is 



9 



3 u 



M° 2 M° 3 =S ab ^gt 



if q a , qb, q c and qa are either all quarks or all anti-quarks, and 



1 2 



$ac g g s p 



M?M2 = 6 ab 6 ae -gi 



2 Q / „",2 £-.1 

•0| 



1 2 



2 //- 



It 



ir 



9 



3 §t 



M 2 M° 3 = -gt5 ab [^-5 ac -j 



ut 



lit 



for any combination of the type q a qb — > g c 9d- 
2) — > processes 

E(5 A , 5b)" = i fc±c) fe±a) f 9 /S B (x b , k ±b ) X 

M°| 2 + |M 2 °| 2 ) + P q z P g z (\M°\ 2 - |M 2 °| 2 )} 



(68) 



(69) 



- A D c/q1 (z, k ±C ) fg/S A (Xa, k± a ) fg/S B (&6, &±b) 



X 



(70) 



M° M 2 ° [P« sin(v9i - v?2 + 0c) " ^ cos(v?i - y? 2 + 0c) 



In this case the amplitude M§ is zero because it violates helicity conservation, and 



\M[ 



0|2 



~ 9 9s 



9s 2 \ 



|2 



- 9 9 S 



u 9 m 2 
"I + II 2 " 



A^A^ = -^ 1-1 + 



9 us 

IF 



(71) 



3) — > gg processes 



2(5,4, 5b) 9 ™ = 2 ^/ g (+ &_i_c) / ff /5 A K, k ±a ) f q/ s B (x b , k ±b ) x 
M°| 2 + |M 3 °| 2 ) + P g z P q z (\M°\ 2 - |M 3 °| 2 )} 



\ AN Dc/ql( z , k ±C') fg/S A (Xa, fc_La) fq/S B ( x b, k ±b) 



X 



M° M 3 ° [P* sin(^ - ^3 + 0c) " ^ C0S(V9! - V?3 + 0c) 



(72) 



where now the amplitude M 2 is zero because of QCD helicity conservation and the ampli- 
tudes [M^gg^gg 8,11(1 [Mg] gg^gg C&I1 be Obtained frOIIl [M?]qg^qg 8,11(1 [M^gg^gg bj appl^g 

Eq. (64). 
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4) qg — > processes 



£(SU, S B ) q9 ^ 9q = - D C / g {z, k± c ) f q /s A (x a , k±a) fg/S B (Xb, k ±b ) X 



1 2 
3 I 



jfjMiT + |m 3 °i 2 ) + p q z p 9 z (iMiT - |m; 

+ \ ^Dc/Tfiz, k ±C ) fq/S A ( x a, k ±a) fg/S B ( X b, k ±b) X 

I M° M 3 ° [T/ cos(v9i - y? 3 + 20g) + T/ sm(<^ - y? 3 + 20^) 



(73) 



where again the amplitude M 2 ° its ze±u ueuaus 
[M^] qg _, gq and [M"^-^ can be obtained from 

5) gq — > processes 



is zero because of helicity conservation and the amplitudes 
obtained from [M^] qg _, qg and [M 2 °] by applying Eq. (65). 



£(£4, 5b) 9 ™ = i £ C/ff (z, fc±c) / fl /5 A (^a, *!±a) / ? /S s K k ±b ) X 

{(iM?) 2 + |M 2 °| 2 ) + P g z P q z (\M°\ 2 - |M 2 °| 2 )} 

+ \ A N D C /T?(Z, k± C ) fg/S A ( x a, k ±a) fq/S B ( x b, k ±b) X 

j^ M 2 ° [T 9 cos(v?i - y? 2 + 20^) + T 9 sin(v?i - y? 2 + 20^)" 



(74) 



where the amplitude M 3 is zero because it violates helicity conservation and the amplitudes 
Wi} gq ^gq and [M$] gq ^ gq can be obtained from [M^] gq ^ qg and [M$] gq ^ qg by applying Eq. (65). 

6) qq -> # c # d processes 

£(SU, = 1 D c/ff (z, fc ±c ) f q/SA (x a , k ±a ) fq/s B (x b , k ±b ) x 



(1 - P* P*) (\M°\ 2 + \M°f) 



(75) 



+ 2M 9 ° M° 



2 -'"a 



^ + P'y PI) C0S(^3 " *> 2 ) " (P« PI - P' Pf) Sin(^ 3 - 



where the amplitude M° is zero because of helicity conservation and 



27 
64 



t 9 f 



64 



i 2 , 



(76) 
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The expression for qq — > gg is obtained from (75) with the replacements q <-> q and y9 2 <-> y?3- 

7) £ a 06 gg processes 

{(l-n a ^)(|M 2 | 2 + |M 3 °| 2 ) (77) 
+ 2M 2 ° M 3 ° [(T^ 7* + T 2 a T 2 b ) cos(y? 3 - </? 2 ) - (l? T b - T 2 a T b ) sin(^ 3 - } , 

and the relevant amplitudes are the same as in Eq. (76), multiplied by the factor 9/64. The 
expression for gg — > qq is obtained from Eq. (77) with the replacements q <-> q and </? 2 <-> y?3. 

8) # a #6 -> processes 



S (^,^)~ = , A,,.(*,^) W^«) W*,*x>) x 

{(|M°| 2 + |M 2 °| 2 + |M 3 °| 2 ) + P 2 a P 2 6 (|M?| 2 - |M 2 °| 2 - |M 3 °| 2 ) 



+ 2M 9 ° M° 



2 ^ 



(7? T+ + T 2 a T 2 b ) cos(^ 3 - + (T 2 a + T« T 2 fe ) sin(y? 3 - ^2) 



+ g A D C/Tf(z, k ±C ) fg/S A (Xa, fc_La) fg/S B (x b , k ±b ) X 

I M° M 2 ° fr^ cos(v?i - y? 2 + 20^) + T 2 a sin(v?i - y? 2 + 20g) 



(78) 



+ M° M 3 ° 



7? cos(v?i - <f3 + 20^) + T 2 b sin(v9i - y? 3 + 20^) 



where 



2\ 



|M 2 T 



|M 3 T 



M 2 °M 3 ° 



9 4 ^ 
2 9s s 2 



u u 

1 + 7 + F 
.'/„■ — I 1 + 4 + ^ 

M M 2 ^ 



r 4 
2"- s 



9 ,t 2 
2 
9 
2 
9 



(79) 
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IV. POLARIZED CROSS SECTION AND SPIN ASYMMETRIES 



Knowing the kernels £(Sa, Sb), we could now proceed with the computation of any polar- 
ized cross section and spin asymmetry, according to our spin and k^ dependent factorization 
scheme, 

Ec da (A,s A )+(B,s B )^c+x 



d 3 p. 



a 



= E 



a,b,c,d 



dx a dx b dz (fk ±a d 2 k ±b d 3 k ±c $(k ±c • p c ) J(k ±c ) 

\K)lX l X a X h Z l S 
\ab^cd 



x Z(S A , S B r^ cd (x a , x b , z, k ±a , k ±b , k ±c ) 5(s + t + u), (80) 



where the sum over all kinds of partons leads to the 8 kernels £(<Sa, Sb) explicitely given in 
Eqs. (67)-(79). 

In the remaining of the paper we shall consider the unpolarized cross section and the 
transverse single spin asymmetry A N and show numerically how much different effects can 
contribute to their values. The single spin asymmetry A at, measured in p^p — > nX scatter- 
ings, is defined as 

d<jl — da^ , . 

A N = , ; , ; , , (81) 



dal + da^ ' 

and requires the evaluation and integration of the quantities E(|,0) — £(|,0) in the nu- 
merator, and £(T,0) + £(|,0) in the denominator. Indeed, the difference and sum of these 
kernels have to be evaluated for each elementary process ab — ■> cd: we shall explicitly show 
the analytical formulae corresponding to four channels only, which serve as examples; all 
the other contributions can be straightforwardly computed in a similar way. 

For the numerator of the single spin asymmetry, for the process A^ B — > C X, we consider 
explicitely the following channels: 

N-a) q a q b -> q c q d 



[E(T,0)-E(|,0)]^ 



\ Af a/A 1 (X a , k±a) fb/B^Xh, k 



+ 2 A~f: 



e s y tf( X a, k± a ) C0S(V9 3 ~ (fi) ~ Af[ 

xAf b Sy/B (x b ,k ±b )M° 2 M° 3 ^ ' ' * 



; ±b ) [|M 1 °| 2 + |M 2 °| 2 + IM3T 



+ 



r lD c /c(z, k ±c ) 
A~f% n {x a , k ±a ) cos(v?i - y? 2 + <t>c) ~ A f) 



} c/c(z, k ±c ) 

l x /l(Xa,k± a ) Sin(^ 3 - <p 2 ) 



(82) 



i> H c) 



p s x /l(x a , k ±a ) sin(v?i - y? 2 + 0c 

x fb/B (x b , k ±b ) M° M 2 ° A N D c/C , (z, k ±c ) 
+ X - Af a/A r(x a , k ±a ) Af b Sy/B (x b , k ±b ) cos( V 9 1 -<p 3 + <j%) M° M 3 ° A N D c/c ,(z, k ±c ) 
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N-b) qq -> g c g d 



[E(t,0) -E(|,0)] 9 ™ 9d = 

- A/ ?m t (s„ fc ± ,) f q/B (x q , k± q ) [ |M 2 °| 2 + |M 3 T] fc±c) 

+ 2 [A - /* B/T (a:„ fc ±g ) cos(y? 3 - ^2) - A fl/^(x q , sin (v?3 - V2) 
x Aj* fe^) M 2 ° Ml b c/g {z, k ±c ) 



(83) 



N-c) qg 



qg 



[E(t,0)-E(|,0)p- 



•9S 



+ 



- Af q/Al (x q , k ±q ) f g/B (x g , k ±g ) [ \M»\ 2 + |M 2 YJ Ay,^, k ±c ) 
A ~flrt( x q, k u) cos (v?i - ¥2 + <f>c) - A fln( x n k ^ sin (^i - V^2 + 0c) 



x k ±g ) M° M 2 ° A N D c / q i (z, k ±c ) 



(84) 



N-d) 5^ -> 5- c 5( d 

[E(t,0)-E(|,0)] 9a " c9d = 

jA/ s/i i(x a) fc lfl )/ 9/Jj (x 6 ,M [|M{f + |M 2 T + |M 3 °| 2 l D c/g (z,k ± c) 



+ 2 



A /|- l/T (x a , fc ±a ) cos(v? 3 - y? 2 ) + A/f- 2/T (a; a , fc ±a ) sin(y9 3 - <p 2 ) 



+ 



x Aft l/B (x b , k ±b ) M 2 ° M 3 ° £> c/ff (z, fc ±c ) 
^~/fi/t(^o,fe±a) cos(</?i -ip 2 + 2(f)%) + Afl,Jx a , k± a ) sin(v?i - v? 2 + 20c) 



/s/sfo, fc ±b ) M° M 2 ° A N D C/Tf (z, k ±c ) 



(85) 



+ - A/ g/AT (x a , k ±a ) Af 9 Ti/B {x b , k ±b ) cos( V 9 1 - v? 3 + 20g) M° M 3 ° A N D c/7 * (z, k ±c ) . 

The above 4 cases have been obtained respectively from the kernels in Eqs. (67), (75), 
(74) and (78), taking into account that 

P x fa/1 = ~ P x fa/ 1 ^2 9 fg/\ = fg/i 

pa v Ln - Py h,i = 2 A-f Sy/] T*f gn - Tff g/l = 2 A-/f /T 

P b x h/B = Af Sx/B = P b z f b /B = Af Sz/B = 7?f g/B = Afr 2/B = 

as one can see from Eqs. (B23), (B46), (B14), (B13) and (B40). 

Let us inspect Eq. (82), which has an immediate partonic interpretation. The first line 
contains the Sivers effect, where the Sivers distribution function for quark a appears in 
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association with the unpolarized parton distribution function (PDF) for quark b, with the 
unpolarized elementary cross section and with the unpolarized fragmentation function (FF) 
for quark c; the second and third lines correspond to the Boer-Mulders effect, in which the 
Boer-Mulders PDF for quark b is convoluted with a complicated combination of distribution 
functions for quark a which, once integrated over the intrinsic transverse momentum kj_ a , 
is somehow related to the transversity function Arg(x a ) or h1(x a ), and the unpolarized FF 
for quark c; the fourth and fifth lines contain the Collins term, coupled to transversity 
distributions, which was already extensively discussed in Ref. [2] (notice that, exploiting 
Eqs. (60), (62) and (B25), one can explicitely show that this term exactly agrees with that 
in Eq. (56) of Ref. [2]); finally the sixth line contains a "mixed" term in which all three 
effects (Sivers <g> Boer-Mulders <g> Collins) appear together. 

Notice that Eq. (85), corresponding to gg — > gg elementary scattering, has the same 
structure as Eq. (82), related to qq — > qq elementary channel: while the Sivers function can 
be defined also for gluons (first line) the other terms correspond to linearly polarized gluons 
inside an unpolarized hadron ("Boer-Mulders-like"), to distributions of linearly polarized 
gluons inside a transversely polarized hadron ( "transversity-like" ) and to the fragmentation 
of linearly polarized gluons into an unpolarized hadron ("Collins-like"). 

In Eqs. (83) and (84), related respectively to qq — > gg and qg — > qg elementary scatter- 
ings, one can recognise the Sivers contribution (first line), the (transversity ® Boer- Mulders) 
[second and third lines of Eq. (83)] and the (transversity (g> Collins) [second and third lines 
of Eq. (84)] effects. 

Concerning the denominator of the single spin asymmetry, da^ + da^ = 2 da unp , the 
relevant quantity we have to calculate is [S(f, 0) + S(|, 0)] for each partonic process ab — > cd. 
We present explicit results for the same channels we have considered above. 

D-a) q a q b -> q c q d 

[£(T,0) + £(|,0)]^™ d = 

fa/A(x a ,k ±a )f b/B (x b ,k ±b ) [ | M° | 2 + | M 2 ° | 2 + | M 3 °| 2 ] D C /c(z,k ±c ) 
+ 2 Af: y/A (x a , k ±a ) Af b Sy/B (x b , k ±b ) cos(v? 3 - V? 2 ) M 2 ° M 3 ° D c/c (z, k ±c ) (86) 
+ [fa/ A {x k ±a ) Af b Sy/B (x b , k ±b ) cos( V 9 1 - v? 3 + <Pc) Ml Ml 



+ Af: y/A (x a , k ±a ) f b/B (x b , k ±b ) cos(v?i - V 2 + 4>c) M^M Q 2 I A JV D c/c1 (z, k ±c ) 
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D-b) qq -> g c g d 



[S(t,0) + S(|,0)] 9 ™ Sd = 

fg/A(x q , k ±q ) f q/B (x q , fcj.,) [ |M°| 2 + |M°| 2 ] A?/^*' fc±c) 



(87) 



+ 2 A/? fe ±9 ) Af /B (x q , k ±q ) cos(y? 3 - y? 2 ) M 2 ° M 3 ° £> c/g (z, fc ±c ) 



D-c) qg 



qg 



[E(t,0) + E(|,0)]" = 

/#(^„ k ±q ) f g/B (x g , k ±g ) [ \M°\ 2 + |M 2 °| 2 ] £c/«(*, fc±c) 
+ Af sy/A (x q ,k ±q )f g/B (x g ,k ±g ) cos(ip 1 - V2 + ^)M 1 M 2 A N D c/q ,(z,k ± c) (88) 

D-d) -> fi^d 



[S(t,0) + S(|,0)] 9a " c9d = 

fg/A(x a , k± a ) fg/B(x b , ^_Lf>) 



■0(2 



|M 2 °| 2 



|M 3 °| 2 



(89) 



+ 2 Afe l/A (x a , k ±a ) Af h Tl/B {x b , k ±b ) cos(v9 3 - V? 2 ) M 2 ° M 3 ° £> C / g (^, fc±c) 

+ [jg/A^a, fc±a) Aj\ /B (x b , k ±b ) C0s(</?i - V? 3 + 20^) M° M 3 ° 

+ Afe l/A (x a , k ±a ) f g/B (x b , k ±b ) cos(v9i - y? 2 + 2 0g) M° M 2 ° ] A N D C/T ?(z, k ±c ) . 

As for the asymmetry numerator, Eqs. (86) and (89), corresponding respectively to the 
elementary scatterings qq — > qq and gg — > gg, have the same overall structure. In this case, 
the first line corresponds to the usual unpolarized term, the second line to a double Boer- 
Mulders (or "Boer-Mulders-like" ) effect, whereas the third and fourth lines contain a mixed 
term in which the Boer-Mulders and Collins (or "Boer-Mulders-like" and "Collins-like") 
effects appear together. Regarding the elementary qq — > gg and qg — > qg channels, Eqs. (87) 
and (88) show that there are two terms contributing to the unpolarized cross section: the 
usual unpolarized term in both cases, the double Boer-Mulders effect for qq — > gg and the 
mixed (Boer-Mulders ® Collins) effect for qg — > qg. 

It might be surprising to notice that several spin and k± dependent mechanisms could 
also contribute to the unpolarized cross section. Their numerical relevance will be studied 
in the next Section. 
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V. NUMERICAL ESTIMATES OF MAXIMAL CONTRIBUTIONS OF SINGLE 
TERMS 



We have now explicit and comprehensive analytical formulae to compute cross sections 
and spin asymmetries, coupling LO QCD interactions and soft physics; all this basic in- 
formation is contained in the kernels, as given in Sections III and IV, to be inserted into 
Eq. (80). 

These kernels, their differences and sums, contain many unknown functions (the soft 
part), which we have interpreted in terms of parton polarizations and distribution or frag- 
mentation amplitudes; we also give their expressions according to the notations of the Am- 
sterdam group (Appendix C). In summary, there are, for each kind of partons, 8 different 
distribution functions and 2 fragmentation functions (into unpolarized hadrons). Out of 
these, only the unpolarized PDF, the helicity distributions and the unpolarized fragmenta- 
tion functions (at least for pions) can be considered as rather well known, from experimen- 
tal information gathered in inclusive and semi-inclusive Deep Inelastic Scattering processes. 
Some approximate information has been very recently extracted also on the quark Sivers 
distribution [15-17] and Collins fragmentation functions [17]. 

This lack of information might induce the thought that any realistic evaluation of physical 
observables, through the scheme of Eq. (80), is hopeless; nevertheless, such a scheme, in 
simplified versions, has already been successfully used to compute transverse single spin 
asymmetries [1, 9] and unpolarized cross sections [1]. Actually, its spin and correlations 
are unique in order to understand and predict many observed and measurable spin effects. 

The way out of this worrying thought is naturally offered by the very structure of our 
scheme and its exact kinematical formulation, and was already partially explored, concerning 
the contribution of the Collins mechanism, in Ref. [2]. The many phases appearing in 
the elementary interactions, due to the noncollinear configurations [see Eq. (60)], once the 
integration over the non observable intrinsic motion is performed, lead to large cancellations 
of most contributions from the new unknown functions. One can realize that looking, for 
example, at Eq. (82) which shows the qq — > qq contributions to the SSA A N . Apart from 
the first line (the Sivers mechanism) all other terms contain the complicated y?i,y?2 and (p 3 
phases, whose integration almost completely cancel their numerical values. 

These effects can be shown in a quantitative way. We have evaluated the different contri- 
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butions to the SSA An and to the unpolarized cross section, taking for each of the unknown 
functions their upper bounds, originating from basic principles (like |P?| < 1). We have 
indeed largely overestimated each single contribution. More precisely, we have adopted the 
following strategy: 

• We have followed Ref. [1], assuming a gaussian k± dependence for all distribution 
functions, with ^ (k\) = 0.8 GeV/c; the same k±c dependence as in Ref. [1] has 
been assumed for the fragmentation functions. At relatively low pt the inclusion of 
k± effects might result in making one or more of the partonic Mandelstam variables 
smaller than a typical hadronic scale. In this case perturbation theory would break 
down. In order to avoid such a problem and extend our approach down to px around 
1-2 GeV/c, we have introduced a regulator mass, fx = 0.8 GeV, shifting all partonic 
Mandelstam variables, that is 

i /i 2 , ii^u- /j, 2 , s -> s + 2/i 2 . (90) 

Concerning the potential ambiguity in the behaviour of the strong coupling constant, 
a s (Q 2 ), in the low Q 2 regime, we adopt the prescription originally proposed by Shirkov 
and Solovtsov [18]. As renormalization and factorization scales Q = p* T j1 is used, 
where p^ is the transverse momentum of the fragmenting parton in the partonic cm. 
frame. A comprehensive study of these and related aspects can be found in Ref. [1]. 

We only stress that, even if the magnitude of each contribution to the unpolarized 
cross sections, in particular at the smallest px values, is sensitive to these choices, 
their relative magnitudes (which are being studied here) are almost not affected at 
all. Moreover, for the SSA (which is a ratio of cross sections) this dependence is 
definitely strongly reduced. 

• The unpolarized PDF have been taken from Ref. [19] and the fragmentation functions 
from Ref. [20]. We have used Eq. (80), taking into account all its partonic contribu- 
tions, and not only those shown as an example in Section IV. 

• All unknown polarized distribution functions have been replaced with the correspond- 
ing unpolarized distributions. In some cases this is certainly an overestimate: for the 
transversity distribution it violates the Soffer bound [21]. 
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• The Sivers and Collins functions have been chosen saturating their positivity bounds: 

A N fa/At (x a , k ±a ) = 2 f a/A (x a , k ±a ) A N D c/q , (z, k ±c ) = 2 b c/q (z, k LC ) ■ (91) 

• Whenever different pieces could combine with different signs (e.g., Sivers or Collins 
functions for different quark flavours), we have summed them assuming the same sign, 
in order to avoid any kind of cancellations not resulting from phase space integrations. 

Our results are shown in Figs. 1-5, and we shortly comment them. 

Fig. 1 shows the different contributions to the unpolarized cross section, see Eqs. (86)- 
(89) for guidance. Our numerical estimate is performed for pp — > it X, in the kinematical 
region of the E704 experiment. Our result clearly proves that the usual contribution involv- 
ing f a /A <S> fb/B® Dq/c largely dominates; even assuming the polarized distributions as large 
as the unpolarized ones and summing additively all of them, their final contributions to the 
unpolarized cross section, after integration over all intrinsic k±, are at least one order of 
magnitude smaller. 

Fig. 2 shows again the different contributions to the unpolarized cross section, for pp — > 
7T° X processes, in the kinematical region of the STAR experiment at RHIC. The dominance 
of the usual f a /A <S> fb/B <8> D c / C term, in comparison with all other contributions, is clear 
again; the second most important contribution, the Boer-Mulders (g> Collins term, is one 
order of magnitude smaller. 

In Fig. 3 we plot the different maximised contributions to A N , for the E704 experimental 
configuration and p^p — > n + X processes, for which very large values of A N have been 
measured [22]. One sees that the Sivers mechanism is largely dominant, that some effects 
might originate from the Collins function and all other contributions are negligible. Notice 
that while the Sivers effect is maximised only in the choice of the Sivers function, the Collins 
contribution is maximised both in the choice of the Collins function and the transversity 
distribution. We have shown separately the quark and gluon Sivers contribution; there 
might be a negative Xp region where one could eventually gain some information on the 
(maximised) gluon Sivers function. 

In Fig. 4 we plot the different maximised contributions to An, for the kinematical region 
of STAR-RHIC experiment, which also has measured non zero values of A^ in p^p — > 7r° X 
processes [23]. Again, the Sivers mechanism gives the largest contribution, some effects 
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"w 10" 35 r Pj=1-5GeV/c 
10" 36 r ^s = 19.4 GeV 
10" 37 r (E704) pp->ji°X 



0.2 0.4 0.6 0.8 1 

X F 

FIG. 1: Different contributions to the unpolarized cross section, plotted as a function of xf, for 
pp — > 7T° X processes and E704 kinematics, as indicated in the plot. The three curves correspond 
to: solid line = usual unpolarized contribution; dashed line = Boer-Mulders £g> Collins; dotted line 
= Boer-Mulders (g> Boer-Mulders. 

might remain from the Collins mechanism and all other contributions are negligible. At 
negative xp all contributions are vanishingly small. 

In Fig. 5 we plot the different maximised contributions to A^, for the kinematical region 
of the proposed PAX experiment at GSI [24], p^p — > n + X . The situation is similar to that 
for the E704 case, with the difference that there might be, at large negative ij?, a region 
where the (maximised) gluon Sivers function gives a sizeable contribution. 



VI. CONCLUSIONS 



We have discussed in great detail a QCD based hard scattering formalism to compute 
unpolarized and polarized inclusive cross sections for the production of large pt particles 
in hadronic interactions. In the absence of rigorous results, we have assumed a factorized 
scheme in which long distance non perturbative physics and short distance pQCD interac- 
tions are separated and convoluted; such a factorization has been proven in collinear QCD, 
but has to be considered as a model when intrinsic motion of partons - effectively introduc- 
ing higher-twist effects - is allowed for. This is the first study in which the intrinsic k± of 
all participating partons is taken into account. This intrinsic motion of partons, generated 
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Vs = 200 GeV 
(STAR) pp->rc°X 



0.2 0.4 0.6 0.8 1 

X F 

FIG. 2: Different contributions to the unpolarized cross section, plotted as a function of xf, for 
pp — > 7T° X processes and STAR kinematics, as indicated in the plot. The 2 lines correspond to: 
solid line = usual unpolarized contribution; dashed line = Boer-Mulders £g> Collins. The Boer- 
Mulders ® Boer-Mulders contribution is not even noticeable at the scale of the figure. 
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°-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

X F 

FIG. 3: Different contributions to An, plotted as a function of xp, for p^p — > ir + X processes 
and E704 kinematics. The different lines correspond to: solid line = quark Sivers mechanism 
alone; dashed line = gluon Sivers mechanism alone; dotted line = transversity 8) Collins. All other 
contributions are much smaller. 
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0.6 



0.4 



0.2 



(STAR) pp->ji u X 
Vs = 200 GeV 
y = 4 
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0.8 



X F 



FIG. 4: Different contributions to An, plotted as a function of xp, for p'p — > tt° X processes 
and STAR kinematics. The different lines correspond to: solid line = quark Sivers mechanism 
alone; dashed line = gluon Sivers mechanism alone; dotted line = transversity (g) Collins. All other 
contributions are much smaller. 
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°-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

X F 

FIG. 5: Different contributions to An, plotted as a function of xp, for p^p — > ir + X processes and 
PAX kinematics, as indicated in the plot. The different lines correspond to: solid line = quark 
Sivers mechanism alone; dashed line = gluon Sivers mechanism alone; dotted line = transversity (g> 
Collins. All other contributions are much smaller. 
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both by confinement and QCD dynamics, plays little or no role in unpolarized processes 
at very large energy, when all relevant momenta are much higher than the average (k±); it 
is however crucial in unpolarized processes at intermediate energies [1] and, even more so, 
in the understanding of spin effects and polarized phenomena. For these, partonic spin-fc^ 
correlations are of fundamental importance: an ever increasing number of spin experiments 
and spin measurements is proving that [22, 23, 25]. 

Eq. (1) is our central point; it is essentially a QCD parton model, in which LO (in 
a s ) pQCD interactions couple to parton distribution and fragmentation functions; intrinsic 
motion is fully taken into account in soft physics and in the elementary interactions. As it 
is well known, this allows new soft partonic functions which would vanish in the collinear 
limit; however, it also introduces in the hard partonic interactions many fcj_-dependent 
phases, which strongly affect the convolution of the soft and hard parts. Luckily, it proves 
that such complicated convolutions involving many phases and many soft functions, have the 
simplifying result of strongly suppressing most contributions to (A, Sa) + (B, Sb) — > C + X 
inclusive processes. Concerning transverse single spin asymmetries An, this leaves at work 
essentially only one spin-fcj_ correlation, namely the Sivers mechanism [11]. This allows to 
explain many measured and intriguing values of A^ [1, 15]. 

We have fully discussed all soft functions, with attention to their physical partonic inter- 
pretation, both in terms of polarized distribution and fragmentation functions and in terms 
of the amplitudes relating partonic and hadronic properties. We have also explicitely shown 
the exact relationships between different notations widely used in the literature; this should 
help in understanding and using the fcj_-dependent factorized scheme. Then, we have numer- 
ically shown the suppression of many contributions, both to the unpolarized cross section 
and the SSA An- This confirms and completes the work of Ref. [2]. 

Many more applications of Eq. (1), modified to hold for different processes, can easily be 
foreseen. This has been done concerning the Sivers asymmetry in SIDIS processes [15] and 
can be extended to the SIDIS Collins asymmetry [26] ; single and double spin asymmetries in 
single particle inclusive production and Drell-Yan processes can equally well be studied, and 
so on. Some information on Sivers and Collins functions is already available from ongoing 
experiments [25, 27] and more is expected; a consistent understanding and computation of 
high energy spin effects, in the framework of a factorized QCD based model, is building up. 
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APPENDIX A: DETAILED k ± KINEMATICS 

We give here, for completeness and the reader's convenience, a detailed account of the 
partonic kinematics with the full inclusion of all transverse momenta, following Refs. [28] 
and [1]. As throughout the paper, we consider the hadronic reaction A B — > C X in the AB 
center of mass frame with A moving along the positive Z cm -axis and we fix the scattering 
plane as the (XZ) cm plane. We neglect all masses, both the hadronic and the partonic ones. 

The 4-momenta of hadrons A, B, C then read 

^ = ^(1,0,0,1) ^ = ^(1,0,0,-1) p£ = (E C , Pt, 0, p L ) , (Al) 

with E c = tJpt + Pl and s = (p A + Pb) 2 - 

For massless partons a, b inside hadrons A, B we introduce light-cone momentum fractions 

%a = Pi/p\ = (P°a + Pl)/(P°A + Pa), x b = Pb/PB = (Pb - Pb)/(P°B ~ Pb) and the transverse 
momenta k± a , k^ b . Their four-momenta then read 

Pa = X a -r-[l + 7= COS0 a , -= Sin0 a , 1 - 

Z y X a S X a y/ S XaV s x a S 

p V^f, , k ±b 2k ±b , 2k ±b k 2 ±b \ 

Pb = x b^T 1 + -5-, ■= COS0 b , -= Sill 0ft, -1 + -=- , (A2) 

2 \ x£s x by /s x by Js xls) 

where k± a>b = |fcj_ a b | and (p a>b are the azimuthal angles of parton a,b three-momenta in the 
hadronic cm. frame. 

The four-momentum of the fragmenting parton c is given in terms of the observed hadron 
momentum p 1 ^, of the light-cone momentum fraction z = Pc/pt and of the transverse 
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momentum k±c of hadron C with respect to parton c light-cone direction. In the hadronic 
cm. frame, we write in general k±c as: 



k± c = k ±c (sm9 k±c cos(p k±c , sm6 k±c sin0 fexc , cos6> fc±c ) , 



(A3) 



and impose the orthogonality condition kj_c ■ p c = via the (^-function 5{k±c • P c )> where 
p c is the unit vector along the direction of motion of parton c. The parton four-momentum, 
p^ = (E c ,p c ), can then be written as 



Pc 



E c = 



E c 



El - k\ c 



(pr - k_L C sin 9 k±c cos <j) k±c , - k± c sin 6 k±c sin (p k±c , p L - k ±c cos 9 k±c ) 



E c + Je% - k\ c 



2z 



(A4) 



and the orthogonality condition kj_ c ■ p c = implies 

_ Je% - k 

d k ± c6(k ±c - Pc) = k^ c dk^ c dO k ^ c d(j) k ^ c — jr 



x 



Vt sin <f% ±c 

W*XO ~ 4 ±c ) + S(4> k±a - (27T - <P° k±c )) 



(A5) 

C ° S ^ = / • U ^^o<- (A6) 

This allows to perform directly the integration over (f> k±c (notice that there are two possible 
solutions to be considered). 

The factor J(z,k±c) entering our basic factorization formula, Eq. (1), is the Jacobian 
factor connecting the parton c to hadron C invariant phase space, defined as 

(A7) 



d 3 p c 1 w , ,d 3 p c 
- = —J(z, k±c)- 



E c :- 



Ec 



which for collinear and massless particles reduces simply to J = 1. With intrinsic motion, 
for massless partons and hadrons: 

( E c + jEl-kl c ) 2 



J(z, k± c ) = J{k± c ) 



(A8) 



"' HE c -kl c ) 

With the expression of parton momenta given in Eqs. (A2) and (A4) one can calculate 



the partonic Mandelstam invariants: 

S = (Pa+Pb) 2 = X a X b S 



„k± a ki_b . k\k\ b 

1 - 2 cos(0 a - (j) b ) + 2 



(A9) 
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where 



i= {Pa-Pc? = -, (A10) 

z 

u = (Pb-Pc) 2 = -, (All) 

z 

/ T + U\ 

s5(s + t + u) = z5[z + — — , (A12) 



E c + JEc- k\ c 
T = -x a ^s . v 2 = — 

2 V E c ~ k±c 



x 



' ^2 \ 2k ±a 

1 + ) \ E c ~ k ±c 7= cos </) a {pT - k ±c sin 9 ' k±c cos <f) k±c ) 



H 7= &_lc sin a sin 6» fc±c sin fc±c - 1 ^ (p £ - k ±c cos fc±o H (A13) 

X a y/S \ X a S J ) 

Er + \ En - fc^ 



= -x h yjs , — 

2JEI - k\ c 



x 



kj_t, \ i ~n To 2/c 



\j Eq — k 2 ±c ^ cos 6 (pr - fc_Lc sin 9 k±c cos fc±c ) 



H 7= sin b sin 6» fc±c sin fc±c +1 — (Pl ~ k_LC cos 6» fc±c ) \ . (A14) 

^feV s \ x b S / ) 



The phase space integrations must obey some constraints, originating from physical re- 
quests. Besides the trivial bounds < x a; b, z < 1, < ai t < 27r and < 9 k±c < it, we 
require that, even including intrinsic transverse momentum effects, a) each parton keeps 
moving along the same direction as its parent hadron, p a ^ ■ Pa{b) > 0, and b) the parton 
energy is not larger than the parent hadron energy, E a ^ < E A ( B y This implies the following 
bounds 

1 (A15) 



k± a (b)/Vs < min x a{b)) ^x a{h) {l - x a{b) ) 

Analogously, for the fragmentation process c — > C + X we require p c • Pc > and 
Ec < -Ec (both fulfilled by Eq. (A4), where we have consistently disregarded the solution 
E c = E c — yj Eg — k\ c /(2z)). The last constraint implies the following bound on k± c , at 
fixed z values: 



k±c/E c < 1 (z < 1/2) k ±c /E c < 2yjz(l-z) (z > 1/2) . (A16) 

By requiring |cos0° ±c | < 1, see Eq. (A6), we have a further constraint on k±c, at fixed 
9 k±c , namely 

p L cos 9 k±c - p T sin 9 k±c < k ±c < p L cos 9 k±c + p T sin 9 k±c . (A17) 
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The partonic helicity amplitudes are computed according to Eqs. (60)-(62); the explicit 
expressions, in terms of the Mandelstam variables, of the relevant combinations of the M° 
amplitudes are given in Section III of the text. The phases (pi are defined in Eq. (62). For 
processes involving only quarks and antiquaries they read: 

<Pi = + 

V?2 = \(Za-Zb-Zc + Zd)+ia-L-<l>c (A18) 

v^s = \ + & - & + &) - L - L + € ■ 

Similarly for processes involving also gluons. 

All terms appearing in the above phases are discussed and can be found in Ref. [2]; we 
report them here for convenience and self-consistency of the paper: 

cos 9 q sin 9j — sin 9 q cos 9j cos(0 g — (f>j) 
008 ^ = sin^ (A19) 

sinQ = ^^-^O. (A20) 

Sill t/gpj 

All angles refer to the overall AB cm. frame. 9j and <f>j (j = a, b, c, d) are respectively 
the polar and azimuthal angles of the partons, while 9 q and <f> q are the polar and azimuthal 
angle of the vector q = p a + p b . Here and in the next equations 9 qp . (0 < 9 qPj < it) denotes 
in general the angle between the two vectors q and p y 
The £j angles (j = a, b, c, d) are given by 

S = ^ + £ , (A21) 



where 



cos # 9 sin 0' - sin 9 q cos 0' cos(0 9 - <f/A 

cos 6 = — J — A22 

J sm&„ 



-sin^sint^) 
J sin 6>„„' 



cos 9' — cos 0'- cos 9„t „/ 

c °sVj = •/■/•/■ A24 

J sin sin O / v 

. , sin^sin^-^.) 

sin 77, = — — • (A25) 

J sin 6* / „' 
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The primed angles (#', 0') are obtained via 



Pi = Pi 



Pi -9 , o 
+ Pi 



q° + \fq 1 \ \[q 1 , 

where i = a,b,c,d and q^ = (q°, q) = p% + P& ■ 

The last angle appearing in Eqs. (A18) is <fi", given by: 

tan 0// = sin^ sin 

c sin 0£ cos(0£. — 0'J cos 0' a — cos 6^ sin Q' a 



(A26) 



(A27) 



Finally, the angle <pQ appearing in the fragmentation amplitudes, Eq. (46), is given, in 
terms of our integration and overall variables, by: 



tan 



=F 



Pt 



Eh 



k 2 \ 
K ±c N 



fe±c - PL cos 6> fcxc ' 
^ p T sin # fcxc 



tan6> 



fc±c ' 



(A28) 



where the =F signs refer, respectively, to the first and second 5-function terms in Eq. (A5). 



APPENDIX B: PARTON POLARIZATIONS AND DISTRIBUTION AMPLI- 
TUDES 

An alternative simple physical interpretation can be given to the distribution functions 
(Pj fa/A,s A ) = ^fs /s A making use of the helicity amplitudes fx a ,x x -,x A , which describe 
the soft process A — > a + X. This is the approach used in Refs. [2, 9]. Since the partonic 
distribution is usually regarded, at LO, as the inclusive cross section for this process, the 
helicity density matrix of parton a inside hadron A with spin Sa and polarization vector 
P A can be written as 

PkX* f*/AS A ( X a,k±a) = E Pif,t' A t ^A a ,A ;A A H' a ,X x ;X' A 

X A ,X' A J A> X A 

- E pf A \ F x a A :k ' ( B1 ) 

X A' X A 

having defined 

Fx a 'y =i fx a a H< x -a' , (B2) 

A A' A A TX A ,Xx A XA A °' X A' X A ' V > 

where the Y] stands for a spin sum and phase space integration over all undetected 

J X A ,Xx A 

remnants of hadron A, considered as a system Xa, and the JF's are the helicity distribution 
amplitudes for the A — > a + X process. 
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Eq. (Bl) relates the helicity density matrix of parton a, see Eq. (7), to the helicity density 
matrix of hadron A, given by 

where P A = {P^ cos 4>s A , Ft sm 0su> Pl) * s hadron A polarization vector and <ps A its az- 
imuthal angle, defined in the helicity reference frame of hadron A. Notice that, in this 
Appendix, we consider the most general case in which the transverse polarization of hadron 
A can be along any direction <p SA in the XY plane, whereas in Section II A and throughout 
the paper the specific choice was made of fixing the transverse polarization of hadron A 
along the F-axis, i.e. \ — Sy, which corresponds to <ps A = it/2. 

The distribution amplitudes T depend on the parton light-cone momentum fraction x a 
and on its intrinsic transverse momentum k± a , with modulus kj_ a and azimuthal angle a : 

F\ a ,\ XA ,\ A (xa, fe±o) = F\ a ,\ XA ,\ A (x a , fc± ) exp[iA A a ] , (B4) 

so that 

Ftjfaa, k ±a ) = k ±a ) exp[*(A A - \' A )<f> a } . (B5) 

F x °\? (x a , k±a) has the same definition as F\ A \ A ( x a, k± a ), Eq. (B2), with T replaced by J 7 , 
and does not depend on phases anymore. 

The parity properties of J~\ a ,\ x -,\ A (x a , k^_ a ) are the usual ones valid for helicity amplitudes 
in the 4> a = plane [4], 

?-K^ Xa ^ a = V (-1)**--^ (-1)^^ F K ,x Xa ,x a , (B6) 
where r\ is an intrinsic parity factor such that rf = 1. These imply: 

p-K-K _ (]\2(s A -s a ) (i\(\ A -K)H* A -K) pK,K /g 7 \ 

-\ A ,-\' A V > \ ) A A' A A ' V ' 

Notice that, for Sa = 1/2, the factor (— l) 2 ^ - **) i s positive if parton a is a quark and 
negative if it is a gluon; consequently, some parity relations are different according to the 
nature of the parton involved. For this reason we shall treat quark and gluon distribution 
functions separately. 

By applying Eqs. (B5) and (B7) one can see that there are six independent F's: 

F+: , F++ , F+: , F-_+ , F++ , F+- . (B8) 
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These are in principle complex quantities, but F£+ and F + + are clearly moduli squared 
(see Eq. B2), whereas F+r and F+ + are purely imaginary for gluons and purely real for 
quarks, as given by Eq. (B7). This leaves us with eight independent real quantities, which 
are directly related to the eight distribution functions defined in Eqs. (13-15) (for quarks) 
and (30-36) (for gluons), as we are going to show. 

1. Quark sector 

Let us consider first quark partons. Inserting Eqs. (B3) and (B5) into Eq. (Bl), and 
exploiting the parity relationships (B7), yields, for a generic hadronic spin state: 

P^ SA fa/A,S A = \(1 + P a z ) fa/A,S A = \ (/++++ + F++) + \ P* - 

+ P$ [ReF++ cos(<fe A - <f> a ) + ImF++ sin(0 5A - <j> a ) 
P a ^ SA L,A,s A = Ul- P a z ) fa/A,S A = \ {F+f + F++) - \ P£ [FIX - 



(B9) 



P4 



ReF++ cos(0 Sa - (p a ) - ImF++ sin(0s A - <j> a ) 



(BIO) 



P T SA fa/A,s A = ^ (PS - iPy a ) fa/A,s A = i ImF+- + P£ ReF+ 



alA,S A 

2 - 



++ 

sm(0s A - (p a ) 



• (Bll) 



(B12) 
(B13) 



'(Ft! + F"+) cos(0 5a -(P a ) -i {F+Z - F-+) si 
By summing and subtracting Eqs. (B9) and (BIO), one finds 

f a /A,s A = (F+t + F --) + 2 Pt sin(0 5A - <j> a ) 

P a z L/a,s a = Pt (FIX ~ F + l) + 2 P£ ReF++ cos(0 5a - <f> a ) , 

while from the real and imaginary parts of Eq. (Bll), 

P a x fa/A,s A = 2 Pt ReF+- + Pt (FX: + cos(0 5a - a ) (B14) 

Py fa/A,s A = -2 ImF++ + Pt (FX! - F~ + ) sinfe - a ) . (B15) 

The two above equations can be written in a compact form (which we shall use later) in 
terms of the parton transverse spin 



P a x =P$ cos (P Sa 



P; = P£sin0 Sa , 



(B16) 
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where Sa is the azimuthal angle of the polarization vector of parton a in its helicity frame. 
By multiplying Eqs. (B14) and (B15) respectively by cos0 Sa and sin0 Sa and summing, one 
obtains: 

PtL/a,s a = -2 ImF+- sin Sa + 2 P* ReF+" cos Sa 



+ P4 



F+_ COs(0 Sa - <j)s A + 0a) + F+ + COs(0 Sa + <P Sa ~ 0a) ■ (B17) 



Moreover, one can show that the azimuthal angle of P a in its helicity frame, Sa , and the 
same angle measured in the hadronic helicity frame, <j)' Sa , are related by 

1 2N 



0s a = 0L - 0a + C 



k 



±a 



(B18) 



so that, up to such corrections, Eq. (B17) can be written as 

PTfa/A,s A = -2ImF+-sin(0; a -0 a ) + 2P i A ReF++cos(0 / Sa -0a) 

+ P* \F+: cos(0^ a - 4> Sa ) + F~+ cos(0' Sa + 5a - 20 a ) 



(B19) 



Eqs. (B12)-(B15) express the quark polarizations in term of the distribution amplitudes 
F's and the hadron polarization. One finds eight non zero independent soft functions: 



fa/A 
fa/A,S T 
Px fa/A,S L 
Px fa/A,S T 
Py fa/A,S L 
Py fa/A,S T 
Pz fa/A,S L 
Pz fa/A,S T 



fa/A,S L = (F£ + F++) 

(P++ + F++) + 2 ImP++ sin(0 SA - a ) 



= 2 ReF 



++ 



(F+r + cos(0 5a - a ) 

P, a / aM = -2 ImF+- 

-2 ImF+- + (f+I - F~+) sin(0 5A - a ) 



++ r — 



= 2ReF++cos(05 A -0a). 



(B20) 
(B21) 
(B22) 
(B23) 
(B24) 
(B25) 
(B26) 
(B27) 



Notice also that P® f a /A = 0. 

If we fix <ps A = 7r/2 as done throughout the paper and adopt the notations of Eqs. (13)- 
(15), the above equations read: 



P a x f a /A,s Y = *L/Sy = A/,. /T = (F+r + F-+) sin0 a 



(B28) 
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Py fa/A,S Y 
Pz fa/A,S Y 
Px fa/A,S z 
Py fa/A,S z 
Pz fa/A,S z 




Af Sy/SY = A/ VT = -2ImF+- + (f+I - cos0 a 

A/ Sz/5y ee A/ Sz/T = 2ReF++sin0 a 

A/*.,*, ee Af Sx/+ = 2 ReF+- 

Ms v ,s z = A/ Sy/+ = A/ Sy/A = -2ImF+- 

Af Sz/Sz ee A/ Sz/+ = - F_ + +) 

/ aM + \ Af a/Sy = + + 2 ImF++ cos a , 



(B29) 
(B30) 
(B31) 
(B32) 



(B33) 



(B34) 



which gives the exact expressions of Eqs. (13)— (15) in terms of helicity distribution ampli- 
tudes. In particular, Eqs. (B32) and (B34) allow to obtain the expressions of the Boer- 
Mulders and Sivers functions respectively [see Eqs. (15), (17) and (21)]: 



2. Gluon sector 

Thanks to the formal analogy between Eqs. (7) and (25) the expressions of the circular 
and linear polarizations of the gluons in terms of the corresponding helicity distribution 
amplitudes are closely analogous to those obtained for quarks in the previous subsection. 
One should only pay attention to the parity properties appropriate for spin 1 gluons and 
remember that the F's are now the helicity distribution amplitudes for the A — > g + X 
process. 

One finds that Eqs. (B9) and (BIO) hold true also for gluons, while Eq. (Bll), due to 
the different parity relationships, changes into: 



& N far/A = -2ImF+; 
A N f a/A , = 4ImF++. 



(B35) 
(B36) 



Notice also that: 




(B37) 



4-' SA f 9 /A,s A = \(V ~ iV)f 9 /A,s A = ReF++ + iP£lmF£ 
- \ Pi [(Ft= + M<t>s A -<f> a )+i (Ft ~ F +-) cos fe - <Pa)} , (B38) 



where F+_ and F + + are now purely imaginary quantities. 
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As a consequence, Eqs. (B12) and (B13) keep describing the distributions of unpolarized 
or longitudinally polarized gluons inside a polarized hadron, while Eqs. (B14) and (B15) 
modify into: 

Vf g /A,s A = 2 ReF++ + P$ Im (f+I + F~f) sin(0 5A - a ) (B39) 
Vf 9 /A,s A = "2 Pi ImF+- - Pi Im (f+J - F++) cos(0 5a - <j> a ) . (B40) 

Eqs. (B20)-(B27) now become: 



fg/A 


= fg/A,S L = + F++) 




(B41) 


fg/A,S T 


= (F++ + F++) + 2 ImF++ sin 


{<PS A - 0a) 


(B42) 


V fg/A,S L 


= T// fl/A = 2ReF+- 




(B43) 


T\ fg/A,S T 


= 2 ReF+~ + Im (F+r + F~+) 


sin(05 A - a ) 


(B44) 


T2 fg/A,S L 


= -2ImF++ 




(B45) 


% 9 fg/A,S T 


= -Im(F+r-F-+) cosfe 


-0a) 


(B46) 


P 9 z fg/A,S L 


= {Kt - f++) 




(B47) 


P 9 Z fg/A,S T 


= 2ReF++cos(0s A -0 a ), 




(B48) 



while, choosing 4> Sa = n/2 and following the notation of Eqs. (30)-(36), we have: 

T? f g/A ,s Y = A/« /5y = A/« /T = 2 ReF+" + Im (f+_~ + F~+) cos a (B49) 

7? / sM ,s y = A/f 2/5y = A/» /T = -Im (F+r - sin a (B50) 

Pi f g/ASY = Af Sz/Sy = Afl n = 2 ReF+_+ sin a (B51) 

V fg/A,s z = Af ri/Sz = A/« /+ = A/£ /j4 = 2 ReF+" (B52) 

^ fg/A,s z = A/« /& = A/« /+ = -2 ImF+- (B53) 

^ fg/A,s z = Af Sz/Sz = Af Sz/+ = (F++ - F++) (B54) 

f 9 /A,s Y = f 9 /A + \ A/ fl/T = (F++ + F_++) + 2 ImF++ cos a . (B55) 

The Sivers function (B36) can exist also for gluons, while the "Boer-Mulders-like" function 
is given by 

Vf a /A = Mr U A = A/£ l/+ = 2 ReF+" . (B56) 
Finally, in analogy to Eq. (B37): 

A-fli/s. = \ [A/|- l/T - A/J/J = Im (F+T + F-+) cos a . (B57) 
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APPENDIX C: RELATIONS BETWEEN DIFFERENT NOTATIONS 



1. Quark distribution functions 

Let us compare our notations with those used in the formalism of the Amsterdam 
group [10] (see also Ref. [29]), which is widely used. In this formalism the main object, 

A A A' 

corresponding to our F x °\? (x a , k± a ), is the <&(x a , k ±a ) correlator 



$(x a ,k ±a ) = i 



hV>+ + /it 77 r *L 9il H 77 fi'iT 7 



M \ " M 



M / M 



+ ^ 



(CI) 



M 

By appropriate Dirac projections $' r ^ = Tr(r$) one can single out the various sectors 
of distribution functions. In particular, Y = (n_) a ^ a /2 projects out the fi sector (i.e. all 
the distribution functions relative to an unpolarized quark), namely the usual unpolarized 
distribution function fi(x a , k± a ) and the Sivers function fyf(x a , k± a ): 

Tr = fi- sin(0 5A - a ) , (C2) 

where M is the proton mass and n± = l/\/2 (1, 0, 0, ±1). 

Similarly, the projection operator T = (n_) a 7 Q 7 5 /2 gives the g\ sector (i.e. the dis- 
tribution functions corresponding to a longitudinally polarized quark), namely the helicity 
distribution function g® L (x a , k± a ) and the number density of longitudinally polarized partons 
a in a transversely polarized hadron A, called gi^(x a , k± a ): 

Tr 7 5 = Pt 9il + P^ cos(0 5a - a ) . (C3) 

Finally, to obtain the hi sector (i.e. the distribution functions relative to a transversely 
polarized quark), we have to apply the projector T = ^ia^nJ)^ T 2 7 : 



A; 2 

cos(<j) SA - <j)'J h + ^ cos(20 a - 5a - 0' S J 



with 



+ ^ cos(0 a - <P'J Pt h{ L - ^ sin(0' Sa - a )^ , (C4) 



hi = h ±T + §f 2 hi T . (C5) 
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The relations between the F x a \f inclusive cross sections and the Amsterdam group dis- 
tribution functions can straightforwardly be derived by comparing Eqs. (C2), (C3) and (C4) 
with Eqs. (B12), (B13) and (B19) respectively, obtaining: 



M 



M 
kia 
M 



fl(x a , k± a 
9lL(x a ,k± a 



k-i-a u j 

M 



hi (x a , k± a 
hi(x a , k ±a 



M 



YlyfiXa, k 



'La 



7++ 
++ 



+ F++ = f a 



/A 



-2 ImF++ 
2 ReF++ 



2 ReF+~ 



2 ImF++ 



= F 



+- 



IF 



+- 



(C6) 
(C7) 
(C8) 
(C9) 

(CIO) 

(Cll) 
(C12) 

(C13) 



Notice that, according to the most general forward behaviour of helicity amplitudes (see, 
e.g., Eq. (4.3.1) on page 79 of Ref. [4]), one should have the minimal requirement: 



F^(x a ,k ±a = 0) ~ (k ±a f-- x « +x '^\ , 



(C14) 



which is explicit in the above equations. The proton mass M is assumed in Eq. (CI) as a 
reasonable scale for the intrinsic motion k±. 

Combining Eqs. (C6)-(C13) with Eqs. (B20)-(B27) one can obtain the relationships 
between the Amsterdam functions and the quark polarizations. Using Eqs. (C7), (CIO), 
(Cll), (C8) and (C9) respectively into Eqs. (B21), (B22), (B24), (B26) and (B27), yields: 

k± a 

fa/A,S T ~ fa/A -s T = Af a/ s T (x a , k ±a ) = -2 sin(0 5A - (j) a ) /j^(x , k ±a ) (C15) 



Pxfa/A,S L = Af Sx/+ (x a ,k ±a ) 
Py fa/A,S L = Py fa/A = Af Sy / A (x a ,k ±a ) 



k 



±a 



M 
k]_a 
M 



hi (x a , k± a ) 



P a z fa/A,s L = Af Sz/+ (x a ,k ±a ) = g 1L (x a , k ±a ) 



(C16) 

(C17) 
(C18) 



PzL/A,S T = Af s js T (x a ,k ±a ) 



k 



±a 



M 



cos(0 5a - (p a ) gi T (x a , k ±a ) , (C19) 



which shows that the functions f^, hf- L , h±, giL and g^ T have a direct physical interpretation 
in terms of corresponding polarized quark distributions. 
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Instead, insertion of Eqs. (C11)-(C13) and (C5) into Eqs. (B23) and (B25) gives 



Px fa/A,S T 



Af Sx /S T (x a ,k± a ) 



h 1T (x a ,k ±a ) + 




hi T (x a , k ±a ) cos(0 5a - <f> a ) 



(C20) 



Py fa/A,S T 



Af Sy/ s T (x a ,k La ) 

— hi(x a , k± a ) + h 1T (x a , k ±a ) sin(0s A - a ) , 



(C21) 



which shows that hyr and h\ T are combinations of quark polarized distributions. 
2. Gluon distribution functions 

In Ref. [14] Mulders and Rodriguez discussed the twist-two transverse momentum de- 
pendent gluon distribution functions for spin-1/2 hadrons. Their notation is different from 
ours, and it is worth mentioning the relations which link the two different formalisms. 

Naming conventions in Ref. [14] are set as follows: G and AG indicate gluon distribution 
functions which are diagonal in the gluon helicities, i.e. correspond to either unpolarized 
(G) or circularly polarized (AG) gluons. H and AH indicate gluon distribution functions 
which correspond to linearly polarized gluons in either unpolarized or polarized hadrons 
respectively. As for the quark distribution functions, a T or L subscript indicates that 
the parent hadron is either transversely or longitudinally polarized, and a _L superscript 
shows an explicit dependence of the distribution function on the gluon intrinsic transverse 
momentum. 

Indeed, eight such functions exist: 

• G is the usual distribution function of unpolarized gluons inside unpolarized hadrons, 
corresponding to f g/A = F++ + Eq. (B41); 

• AGl is the distribution function of circularly polarized gluons inside a longitudinally 
polarized hadron A, corresponding to Af 9 s ^ + = — Eqs. (B47) and (B54); 

• Gt is the distribution function of unpolarized gluons inside a transversely polar- 
ized hadron, i.e. the gluon Sivers function, corresponding to A N f g / A i; = 4ImF^j", 



Eq. (B36); 
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• AGt is the distribution function of circularly polarized gluons inside a transversely 
polarized hadron, corresponding to A/* , St = 2ReF_jtj~, Eqs. (B48) and (B51); 

• H 1 - is the distribution function of linearly polarized gluons in unpolarized hadrons, 
which corresponds to A/£ /A = ReF+~, Eqs. (B43) and (B52); 

• Hj; is the distribution function of linearly polarized gluons in longitudinally polarized 
hadrons, which corresponds to Af^ 2 j + = im -^++ > Eqs. (B45) and (B53); 

• AHt and AH^ are related to the distribution function of linearly polarized gluons 
in transversely polarized hadrons, AH' T = AH T — (k\ g /2M 2 )AH^. In this case, it 
is difficult to find a precise relation between the two formalisms, but we can say that 
AH T and AH^ play the same role as F^Z and Fj+ , similarly to the quark case [see 
Eqs. (B28), (B29), (C20) and (C21)]. 

Notice that Eq. (C14) is valid for gluons as well as for quarks. 



APPENDIX D: HELICITY FRAMES 

Our physical observables are computed in the AB cm. frame (overall hadronic frame) 
with axes denoted by X cm , Y cm , Z cm . The helicity frame of a particle with momentum p 
along the direction p = (sin 6 cos <p, sin 9 sin ip, cos 6) - as defined in the hadronic frame - 
can be reached by performing on the overall frame the rotations [4] 

R(<p,e,o) = RY>(e)R Zm (<p). (di) 

The first is a rotation by an angle tp around the Z cm -axis and the second is a rotation by an 
angle 9 around the new (that is, obtained after the first rotation) F'-axis. 

This results in the helicity frames with axes along the following directions (expressed in 
the hadronic frame): 

Xa — X cm Y a — Ycm Z a = Z cm (D2) 

for a hadron A moving along +Z cm , 

Xb = X cm Y B = —Y cm Z B = —Z cm (D3) 
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for a hadron B moving along — Z cm , 

~ - Z cm ^ P ^ t ~ ~ /t-a^\ 

x = yxz y = = Z cm x k ± z=p (D4) 

for a generic particle p. Notice that is the unit transverse component - with respect to 
the Z cm -direction - of p, and that it lies in the (xz) plane. 
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